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ABSTRACT 

We present the results of detailed monitoring of the radio emission from the 
Type Ic supernova SN 19941 from 3 days after optical discovery on 1994 March 31 
until eight years later at age 2927 days on 2002 April 05. The data were mainly 
obtained using the Very Large Array at the five wavelengths of AA1.3, 2.0, 3.6, 
6.2, and 21 cm and from the Cambridge 5 km Ryle Telescope at A2.0 cm. Two 
additional measurements were obtained at millimeter wavelengths. This data set 
represents the most complete, multifrequency radio observations ever obtained 
for a Type Ib/c supernova. The radio emission evolves regularly in both time and 
frequency and is well described by established SN emission/absorption models. It 
is the first radio supernova with sufficient data to show that it is clearly dominated 
by the effects of synchrotron self-absorption at early times. 

Subject headings: galaxies: individual (NGC 5194 [M 51]) - radio continuum: stars - 
stars: mass-loss - supernovae: general - supernovae: individual (SN 19941) 
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Introduction 



1.1. Background 



By 1941 iMinkowskil (ll94lf ) had recognized at least two classes of SNe - Type I with 
no hydrogen in their optical spectra and Type II with hydrogen lines in their optical 
spectra. By 1980 radio observations of supern ovae began to play a role with the clear 



detection of SN 1979C flWeiler & Sramek 



radio supernova (RSN), SN 1970G, by 



1980) (there had been only one earlier detected 



Gottesman et al. 



1972 



Goss et al. 



1973 but the 



Westerbork telescope was unable to cleanly separate its radio emission from a nearby HII 
region) and the monitoring of its rad io light curves at multiple radio frequencies for many 

200ol). The radio detectio n and monitoring of SN 1980K 



years thereafter flMontes et al 



quickly followed (IWeiler et al 



1986; 



Montes et al. 



19981 ). However, both objects were 



optically identified as Type II SNe. At that time, little was known about the radio emission 
from SNe and the stati stics with only two objects were very poor. However, the radio 



detection of SN 1983N (Srame k. Panagia. fe Weiler 



authors), SN 1984L (jPanagia et al. 



19841 . a lso labelled SN 1983.5 1 by those 



1986|), and SN 1990B (Ivan Dvk et al. Ill993 ) 



Panagia et al. 



appea red 



(1986) 



to establish radio detections of Type I SNe. Around that time, 
proposed that, from their optical spectra, Type I SNe are not a homogeneous class but split 
into at least two sub-classes. The appearance of SN 19941, the subject of this paper, finally 
yielded a bright, nearby object for detailed study of a "radio loud" Type I supernova. 

Currently, our understanding of the classes of SNe is such that the category of Type 
I supernovae is clearly bifurcate d. One now speaks o f Type la SNe (which still have never 



been detected in the radio, see 



Panagia et al. 



20061 ) as originating in the deflagration 



of small stars, such as accreting or merging white dwarfs or black holes, and Type Ib/c 
and Type II (of which there are also numerous subclasses) SNe as originating in the core 
collapse of massive stars with a ZAMS masses of > 8 M . 
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Particularly si nce the discovery of an association of SN 1998bw with GRB 980425 



( I Sadler et al 



19981 ). the study of Type Ib/c SNe has greatly increased (as explained below, 
the optical spectroscopic differences between Type lb and Type Ic SNe are so slight that we 
treat them as a single class). However, further wor k has found that most Type Ib/c SNe 



(>97 %) are not associated with detectable GRBs ( iBerger et al. 
20061 ) . 



2003a; 



Soderberg et al. 



relatively poorly sampled in the radio (see, e.g., 



2004 



2005 



c SNe since S A 


J 1998 


3W, most are still 


Berber et al. 


2003b 




Soderberg et al. 



20101 ) and, even at the present time, no Type Ib/c SN has been as thoroughly 



studied in the radio as SN 19941 presented in this paper. 



1.2. SN 19941 



SN 1 9941 in M 51 (NGC 5194) wa s discovered at magnitude V = 13 m 5 on 1994 April 



2.17 UT flPuckett & Armstrong 



following quickly (see IAUC 5961). 



1994 ) with independent discoveri es by other observers 



Richmond fc Filippenko 



1994) confirmed the disco very 



Armus fc Mazzarella 



(1994) 



as a supernova (SN) and noted it as being "very blue/ 
obtained moderate-dispersion spectra and noted "broad undulations superimposed on a 
generally featureless continuum." 



Initial confusion arose as to the type identification of the SN with 



Schmidt fc Kirshner 



(jl994j ) first (incorrectly) identifying it as a Type II SN based on unreduced spectra taken 
by J. Peters and G. Bernstein during 1994 April 3.37-3.52 UT. The spectra, however, 
did suggest the SN location in NGC 5194 (M51) by ide ntifying strong Na D lines at the 



recess ion velocity of the galaxy. By 1994 April 4.4 UT, 



Filippenko. Matheson. &: Barth 



(119941 ) suggested from new spectra that SN 19941 was a Type lb SN. However, they noted 
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the "possible presence of weak Hq," and suggested that, although the presence of hydrogen 
would "technically make it a Type lib," the steep drop of the continuum toward the blue 
and near-ultraviolet was more consistent with a Type lb i dentifi cation, as were other aspects 



of the spectrum. Finally, 
was not a Type II SN and 



Schmidt. Challis. fc Kirshner 



Phillips 



(119941 ) confirmed that SN 19941 



( 119941 ) pointed out the spec tral similarities to the Typ e 



Ic SN 1983V. The Type Ic identification was also confirmed by IClocchiatti. et al. 



(ll994T ). 



(Note: Because the spectral differences between the SN Type lb and Type Ic optical classes 
are slight - Type lb show strong He I absorption while Type Ic show weaker He I absorption 
- and ther e are no obyious r adio differences, we shall hereafter refer to Type Ib/c SNe as 



one class.) 



Richmond 



(119941 ) discovered a red star near the SN position on HST PC images 



from 1992 which he suggeste d as a possible progenitor to the SN, although his suggestion is 



disputed by 



Kirshner 



(Il994f ). 



Already on 1994 April 3.17 U 



Rupen et al. 



(jl994J) detected radio emission at the 



position of SN 19941 with the VLA|±| at both A1.3 cm (22.4 GHz) and A3.6 cm (8.4 GHz). 
The radio observations established an accurate position of RA(J2000) = 13 ft 29 m 54^12, 
Dec(J2000) = +47°ll / 30 / . / 4 with an uncertainty of ± O'/l in e ach coordinate with a n offs et 



Morrison fc Argyle 



(119941 ) 



from the radio nucleus of M 51 of 14'/3 east and 12'.'3 south, 
report a slightly different optical astrometric position of RA(J2000) = 13 /l 29 m 54 1 !072, 
Dec(J2000) = +47°ll'30 / . / 50 with uncertainty of ± 0"08 in each coordinate, a difference of 
~ 0'.'5 in RA and ~ fX'l in Dec. This optical/radio position difference might imply that there 
is a significant offset between the centroids of the optical and radio emission (unlikely since, 
at the distance of M51 at 8.9 Mpc, (X'5 is 21.5 parsecs) or, more likely, that there is an offset 
between the radio reference frame and the FK5 based optical reference frame in that part of 



x The VLA telescope of the National Radio Astronomy Observatory is operated by Associ- 
ated Universities, Inc. under a cooperative agreement with the National Science Foundation. 
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the sky. iTully I ( 119881 ) determ ined a distan ce of 7.7 Mpc (assuming H = 75 km s" 1 Mpc" 1 ] 



to NGC 5194. Accepting the 



Tully 



(ll988[ ) method and using H = 65 km s 1 Mpc 1 , we 
will adopt 8.9 Mpc as the distance to SN 19941. Other dist ances to M51 range from 9.6 Mpc 



(Hq = 57 : derived from HIT regions and brighte st stars by 



8.4 Mpc (iFeldmeier. Ciardullo. &: Jacoby 



Sandage fc Tammann 



197J) to 



19971 1 from planetary nebula luminosity. 



Ult raviolet measurements were reported from 



5.13 by 



Rodriguez- Pascual. Monier. fc Wamsteker 



reported on 1994 April 23.2 UT bv iRudv fc Dotan 



UE on 1994 April 3.35, 4.9, and 



994) and an infrared detection was 



(|l994f l in the J, H, and 2.00-2.32 /im 



windows. An unabsorbed soft X-ray luminosity 



upper limit of < 1.5 x 10 38 erg s 1 fo r 



Lewin. Zimmermann. fc Pietsch 



(11994h 



a distance to M51 of 7.7 Mpc was reported by 
from the ROSAT satellite for 1994 May 22 and subsequent X-ray studies yielded the 
detectio n of SN 19941 with ROSAT at 82 days after explosion with (1.82 ±4.2) x 10 38 



erg s 1 ( jlmmler et al. 



19981 ). The CHANDRA X-ray Facility at 2271 and 2639 days after 



explosion reported (1.4 ± 0.7) x 10 37 erg s 1 and (1.2 ± 0.5) x 10 37 erg s \ respectively 



( jlmmler. Wilson, fc Terashima 



20021). 



Immler. Wilson, fc Terashima 



(120021 ) determined a 



best-fit X-ray rate of decline proportional to t~ l and a circumstellar medium (CSM) density 
profile of Pcsm oc r _1 - 9±ai , consistent a constant mass- loss rate/constant wind velocity for 
the SN progenitor which would yield a CSM density profile of Pcsm oc r ~ 2 - 



Prediscovery limits are also available. iRichmond. et al. I ( 119941 ) reports that SN 19941 



was not present to a limiting R magnitude of 



measurement with the VLA by 



Tongue, et al. 



6.2 on 1994 March 28 UT and a prediscovery 



(119941 ) found a 3cr upper limit for A21 cm 



emission from that position of <2.3 mJy on 1994 February 12. 

Due to its proximity (8.9 Mpc; Tully 1988 recalculated for H = 65 km s" 1 Mpc" 1 ) and 
the fact that such an early radio detection was obtained, we quickly started observations 
at all available wavelengths with the VLA. Regular monitoring of SN 19941 with the VLA 
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at AA1.3 cm (22.5 GHz), 2 cm (14.9 GHz), 3.6 cm (8.4 GHz), 6 cm (4.9 GHz), and 20 cm 
(1.4 GHz) then continu ed until the end of 1994 with sporadic measurements thereafter. 



Pooley fc Green 



(119941 ) also monitored the SN extensively at 2 cm (15.2 GHz) with the 
Ryle Telescope. These data are listed in Table [TJ Two additional observations were taken 
with t he Caltech OVRO millimeter interfe rometer at 3 mm (99 GHz) and at 1.4 mm (218 



GHz) (IPhillips. Chandler, fc Rupen 



199J). These are listed in Tabled. 



2. Radio Observations 
2.1. Flux density measurements 

The VLA radio observations reported here were made between 1994 April 03 and 2002 



April 05. VLA phase and flux density ca 



procedures (see, e.g., 



Weiler et al. 



1986 



ibration and data reduction followed standard 



19901 ) . using 3C 286 as the primary flux density 



calibrator and J1418+546 as the secondary flux density and phase calibrator. As noted in 
the footnotes to Table [31 two measurements on 1994 April 03 used the calibrator J1415+463 
and one measurement on 2002 April 05 used the calibrator J1327+434. Table [3] lists the 
measured or interpolated (or extrapolated) flux densities of the secondary calibrators. For 
the phase calibration, J1327+434, J1415+463, and J1418+546 have defined positions as 
listed in the NRAO VLA calibrator catalogue. 

The flux density measurements for SN 19941 are listed in Tables [1] and [2] and plotted in 
Figure [T] with the filled squares representing actual detections and inverted open triangles 
representing 3er upper limits. Also shown in Figure Q] are the best fit model curves discussed 
below in Section [3] and Section 0] with the set of parameters reported in Table 4. 

The flux density measurement error given for the VLA measurements in Table [1] is a 
combination of the rms map error, which measures the contribution of small unresolved 
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fluctuations in the background emission and random map fluctuations due to receiver noise, 
and a basic fractional error e included to account for the normal inaccuracy of VLA flux 
density calibration and possible deviations of the primary calibrator from an absolute flux 
density scale. These final errors (erf) are taken as 

a) = (eS f + al (1) 

where Sq is the measured flux density, o"o is the map rms for each observation, and e = 0.10 
for 20 cm, e = 0.05 for 6 and 3.6 cm, e = 0.075 for 2 cm, and e = 0.10 for 1.2 cm. For all of 
the Cambridge measurements at 2 cm, e = 0.10 was assumed. 



2.2. Spectral indices 

With such extensive data, we can also calculate the spectral index evolution for 
SN 19941 and these are shown in Figure [2] for four pairs of bands (1.2-2 cm, 2-3.6 cm, 
3.6-6 cm, and 6-20 cm). The parameters listed in Table 4 from the best fit light curves allow 
calculation of the model spectral index evolution curves and these are also shown in Figure 
[2J At the earliest times, the observed spectral indices exceed the maximum value spectral 
index (a = 2.5) possible for a purely synchrotron self- absorption source, suggesting that at 
least some thermal absorption is present in SN 19941. 

The overall model fits to the spectral index evolution are quite good, with the only 
exception being the model 1.2-2 cm spectral index at early times, which is systematically 
higher tha n the observed valu es. This deviation, which was found also for the case of 



SN1993J flWeiler et al. 



20071 ). suggests that, at the highest frequencies, the intrinsic radio 
spectrum has a somewhat flatter spectral index than the best fit value determined by the 
longer wavelength data. 
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2.3. Brightness temperature 

This extensive data set can also be used to estimate the intrinsic brightness temperature 
of the SN 19941 radio emission which can be derived from the ratio of the observed radio 
flux density to the solid angle fl of the emitting region (T B oc S b s /Q). 

Since SN 19941 was not radio bright enough, it was not possible to directly measure 
its a ngular size by means o f VLBI observations, as had been possible for SN 1993J (see, 



e.g 



Marcaide et al. 



20091 . and references therein), but the physical size of the SN shock 
can be estimated from the optically measured expansion velocity and distance to the 
SN. For SN 19941, published expansion velocities at early times range between ~16,500 



km s 1 (jFilippenko et al. 



1995 



Millard et al. 



I A 10830, through -20,000 km s" 1 ([Millard et al. 



30,000 km s 1 that is the extent measured by 



1999) from the absorption troughs of He 



1999) as seen for S i II, up to almost 



Filippenko et al. 



(119951 ) for the blue wing 



of the He I A 10830 absorpt ion line. We shall adopt a com p romis e expansion velocity of 



20,000 km s 1 and, following 



Immler. Wilson. &: Terashima 



( 120021 ) . we shall assume that it 



is constant with time. Therefore, the radius r of the emitting region can be estimated as: 



r = 1.73 x 10 14 (t/1 day) cm 



(2) 



which, for an adopted distance D to SN 19941 of 8.9 Mpc, gives an angular size of 



6 = 1.29 (t/1 day) /iarcsec 



(3) 



and an estimated brightness temperature of 



T B = 3.5 x 10 n S , obs (mJy) A (t/1 day)" 2 K 



(4) 
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The resulting brightness temperature evolution is shown in Figure [3] (left). The 



corres ponding brightness temperature evolution curves for SN 1993J from 



Weiler et al. 



( 120071 ) are also shown in Figure [3] (right). 



From these brightness temperature estimates, we can constrain the value of external 
absorption du e to a presupernova stellar wind by requiring that T R not exceed the limit of 



3 x 10 11 K flKellermann fc Paulinv-Toth 



1969 



Readhead 



19941 ) where the radio emission 



is quenched by inverse- Compton scattering. This implies that any value of K2 (see Section 
13. ip higher than ~ 30 would yield Tb curves in strong disagreement with such a physical 
requirement and puts an upper limit on the amount of external thermal absorption by a 
presupernova mass-loss in a dense stellar wind. 



3. Radio supernova models 



3.1. Parameterized model 



Weiler et al. 



( 119861 ) discuss the common properties of radio SNe (RSNe), including 
non-thermal synchrotron emission with high brightness temperature, turn-on delay at 
longer wavelengths, power-law decline with index (3 after maximu m, and spectra l index 



a (S oc v +a ) asymptotically decreasing to an op tically thin value 



by 



1990T) have shown that the "mini-shell" model of 



Weiler et al. 



Chevalier 



(11982a 



Weiler et al. 



(1986 



b), with modifications 



()1990l ). adequately describes previously known RSNe. In this model, the 
relativistic electrons and enhanced magnetic fields necessary for synchrotron emission are 
generated by the SN shock interacting with a relatively high-density circumstellar envelope 
which has been ioni zed and heated by the initial UV/X-ray flash. Figure H] (taken from 



Weiler et al. 



(120021 )) shows a, not to scale, greatly simplified cartoon of the vicinity of a 



supernova explosion. This dense cocoon is presumed to have been established by a constant 



- 12 - 



mass-loss (M) rate, constant velocity (w) wind (i. e., pcsm oc r 2 ) from a massive ( > 8 M G 



0; 



can also be the source of initial thermal 



Chevalier 



(119981 ) has suggested that 



stellar progenitor or companion. This ionized CS1V 
absorption (Free-Free Absorption, FFA) although 
non-thermal, synchrotron self-absorption (SSA) may be a significant source of absorption 
in the early phases of some RSNe. For the case of FFA, the rapid rise in radio flux density 
results from the shock overtaking progressively more of the wind matter, leaving less of it 
along the line-of-sight to the observer to absorb the more slowly decreasing synchrotron 
emission from the shock region. In the case of SSA, the n on-thermal absorption decreas es 



as the emitting region expands and decreases in density. 



Montes. Weiler. and Panagia 



( 119971 ) also allow for the possibility of intervening HII along the line-of-sight to the RSN, 



providing a constant thermal absorption of the radio emission. 



To parameteri ze the basic characte ristic s of radio supernoy a e, we 



RSN discussion of 



Weiler et al. 



((2002) and 



Sramek fc Weiler 



(12003 



adopt the most recent 



and illustrate the 



regions of the super nova being parameterized in a cartoon shown in Figure H] from, e.g. 



Weiler et al. 



(12002) : 



S(mJy) = K x 



« ft-tc 



5 GHz/ Viday 



1-e 



clumps 



r CSM clun 



g ^internal 

^internal 



(5) 



with 



''"external ^CSM homogcncous "^distant; 



(6) 



2 Even though this model has been previously described in papers such as the references 
given above, we list it again for completeness, for the reader's convenience, and for explana- 
tion of the parameters listed in Table H] 
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where 



7csM homog _ K 2 ( g GH J ^ (7) 



-2.1 

^distant = K 4 ( - ^ ) (8) 



and 



^CSM clumps - ^3 (^) (^) 



5' 



(9) 



with Ki, K 2 , K 3 , and K4 determined from fits to the data and corresponding, formally, 
to the flux density {Ki), homogeneous (K 2 , K4), and clumpy or filamentary (K 3 ) FFA at 
5 GHz one day after the explosion date to- The terms i~csM homogcncous and TcsM clumps describe 
the attenuation due to local, homogeneous FFA CSM and clumpy or filamentary FFA 
CSM, respectively, that are near enough to the SN progenitor that they are altered by the 
rapidly expanding SN blastwave. The TcsM homogcncous FFA is produced by an ionized medium 
that completely covers the emitting source ("homogeneous external absorption"), and the 
(1 — e _TCSMclum P s )r c 7 s 1 M j term describes the attenuation produce d by an inhomogeneous 
FFA medium ("clumpy absorption"; see 



Natta fc Panagia 



1984. for a more detailed 



discussion of attenuation in inhomogeneous media). The r^ant term describes the 
attenuation produced by a homogeneous FFA medium which completely covers the source 
but is so far from the SN progenitor that it is not affected by the expanding SN blastwave 
and is consequently constant in time. All external and clumpy absorbing media are assumed 
to be purely thermal, singly ionized gas which absorbs via FFA with frequency dependence 
z/ -21 in the radio. The parameters 5 and 5' describe the time dependence of the optical 
depths for the local homogeneous and clumpy or filamentary media, respectively. 
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Since it is physically realistic and is needed in some RSNe where radio observations 
have been obtained at early times and high frequencies, Equation also includes the 
possibility for an internal absorption terno This internal absorption (r intcrna i) term may 
consist of two parts - synchrotron self-absorption (SSA; r interna i SSA ) , and thermal FFA 
(Tintcrnai FFA ) due to ionized gas mixed with non-thermal emission. 

^internal TinternalgsA ^"internalppA (10) 



/ v \ a ~ 2 - 5 ft-t \ 5 " 

Wna lssA - K, {^-J ^ — I (11) 



WnaW " K 6 ( 5 GH J ^ ^ J (12) 

with K$ corresponding, formally, to the internal, non- thermal (z/*~ 2 - 5 ) SSA and K 6 , 
corresponding, formally, to the internal thermal (v~ 2A ) FFA mixed with non-thermal 
emission, at 5 GHz one day after the explosion date to- The parameters 5" and 5"' describe 
the time dependence of the optical depths for the SSA and FFA internal absorption 
components, respectively. 

Application of this basic parameterization has been shown to be effective in describing 



3 Note that for simplicity an internal absorber attenuation of the form 



l- e - TCSM i: 



r GSM internal 



which is appropri ate for a pla ne^arallel geometry, is used instead of the more complicated 



expression (e.g. 



Osterbrock 1119741 ) valid for the spherical case. The assumption does not 



affect the quality of the analysis because, to within 5% accuracy, the optical depth obtained 
with the spherical case formula is simply three-fourths of that obtained with the plane- 
parallel slab formula. 
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the physical characteristics of the presupernova system, its CSM, and its final stages 
of evolution before explosion for objects ran ging from the two d ecades of monitoring 



the complex radio emissio n from SN 1979C (IMontes et al 



SN 1 998bw (GRB9804251 flWeiler et al. 



SNe flWeiler et al. 



2002 



20001 ). through the unusual 



200 1_b|), to most r ecent 7-ray bursters/Type Ib/c 



Weiler. Panagia. fc Montes 



2003|) 



4. Discussion 

4.1. Application of models to SN 19941 

As was clearly demonstrated by Weiler et al (2007) for the case of SN 1993J, in order 
to obtain a meaningful representation of radio observations of a supernova, a model must 
be able to account simultaneously for three different aspects, namely: (a) the light curves 
measured at several epochs and in several distinct radio bands, (b) the corresponding 
spectral indices computed for pairs of bands at adjacent frequencies, and (c) the brightness 
temperatures, which at all times should be lower than the well established limit of 3 x 10 11 K 
(Kellermann & Pauliny-Toth 1969, Readhead 1994). 



4.2. SN 19941 radio light curves 

We have excellent time coverage of the radio emission from SN 19941 at 5 frequencies, 
from 22 GHz down to 1.4 GHz. That enables us to make a detailed study of its radio 
properties. In particular, the rising branch of the light curve at early times, where 
absorption (r >> 1) is still important, is well defined in all five main bands. Also, the flux 
density evolution is well characterized by a power law decay at late times, in the optically 
thin phases. Therefore, we are in an excellent position to identify the main processes 
operating in this source and determine their parameters. Actually, the best fit curves as 
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described by the parameters in Table H] and shown in Figure [T] agree quite well with the 
observations, with an overall \ 2 = 4.42 per degree of freedom. 

From this rather good parameterization of the radio emission from SN 19941, we 
conclude that the radio light curves are not well described by a purely thermal absorption 
model such as that given by Equations to [9J In particular, the requirement for non-zero 
K§ and Kq parameter values implies the existence of significant synchrotron self-absorption 
(SSA). In fact, the early optical depth is dominated almost entirely by the SSA (K 5 ) 
component and the FFA component (K&) makes a significant contribution only as the 
emission bec omes more optica lly thin and the brightness temperature begins to drop. As 



predicted by 



Chevalier 



(119981 ) there does not appear to be a highly-structured, thick-CSM 
present since there is no need for a large K 2 or non-zero K 3 term. The necessity of a 
K4 term (i.e., thermal absorption by foreground HII gas; see Section I4T51) is consistent 
with recent deep 20 cm observations that indicate diffuse radio emission from the host 



galaxy in the region near the supernova (IMaddox et al. 



20071 ). Further evidence for 



substantial amounts of HII gas is pr ovided by the 



the immediate vicinity of SN 19941 (IMaddox 



2006 



a emission detected on HST images of 



Maddox et al. 



20071 ). 



It is interesting to note that SN 19941 shows a dip around day 27 in its flux density 
at centimeter wavelengths (frequencies from 4.9 GHz to 22 GHz). While less pronounced 
than that noted for the unusual Type Ib/c SN 1998bw, wh ich is related to the 7 -ray 



burster GRB 980425, it occurs on a similar time scale (see 



Kulkarni et al. 



1998 



and 



http://www.narrabri.atnf.csiro.au/public/grb/grb980425/). This early "dip" followed by 
"re-brightening" for SN 1998bw has variously been explained as due to a slowing and 



re-energization of the SN blast wave (ILi and Chevalier 



in the circumstellar absorbing material ( iWeiler et al 



2001a 



1999) or, more likely, as structure 



b) 



Considering the light curves shown in Figure [T] in more detail, even though the fitted 
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model provides a good representation the agreement is never quite perfect, mostly because 
the model is rather simplified and also because there may be systematic deviations from 
our basic assumptions. In order to analyze the residuals between the model and the 
observations, we express them as logarithms of the ratios of observed to model fluxes 
and plot them in Figure [5] (left). A similar analysis helped us in othe r RSNe to reyeal 



unsuspected features suc h as a periodic flux modulation for SN 1979C ( jWeiler. et al. 



1992 



Schwarz &: Pringle 



19961 ) possibly du e to the presence of a companion, as well as eviden ce 



for in teracting winds in SN 1998bw (I Weiler et al. 



20041 ). 



2001bh and SN 200 lig flRvder et al. 



The residuals for SN 19941 are displayed separately for individual bands in Figure [5] 
(left), and merged in a single plot in Figure [5] (right). As expected, they show scatter 
inherent to non-negligible measurement errors and to a simplified description of a 
complicated phenomenon. However, a systematic feature is apparent, essentially with the 
same timing and comparable amplitudes, at all bands. In particular, in Figure |5] (right) 
we can see a somewhat broad trough, with an apparent local minimum, which is about 0.2 
deep in the logarithm and occurring ~ 27 days after explosion. This trough extends from 
~ 10 up to ~ 55 days after explosion. Since this local minimum is present at all frequencies 
with similar amplitude, it must represent a genuine decrease in the intrinsic emission rather 
than being the result of additional, temporary absorption along the line of sight. In a 
spherically symmetric model it would correspond to a temporary decrease of the average 
density of the CSM by ~ 0.1 in the logarithm of the ratio, or a factor of ~ 0.8 decrease in 
the CSM density, with a duration of ~ 45 days. 

Although a temporary decrease in the mass-loss rate from the presupernova star 
cannot be ruled out, we feel it worthwhile to also consider a model in which the SN 19941 
progenitor was a Wolf-Rayet (WR) star in a wide binary system whose companion was also 
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an early type star (most likely a B-type star on the main sequence). Such a binary system 
could create no n-spherical CSM density structures by the interacting winds of the system 



stars (see, e.g. 



Schwarz fc Pringle 



19961 ). If such were the case, for SN 19941 the observed 



minimum in the residuals would correspond to a s eparation between the two stars of 



4.7 x 10 15 cm or ~ 310 AU [see also SN 1979C f lWeiler. et al. 



1992 



Schwarz fc Pringle 



19961 ) where the radio emission modulation is possibly due to the presence o f a companion 



as well as evidence for interacting winds in SN 1998bw (IWeiler et al 



2001b|) and SN 200 lig 



(IRyder et al. 



200J)]. 



4.3. SN 19941 spectral index evolution 

Examination of the spectral index evolution in Figure |2] shows three things: 1) The 
model listed in Table H] generally describes the spectral index evolution quite well except 
at at 1.2-2 cm. Apparently, at higher frequencies the model is beginning to fail, 2) The 
spectral index evolution, particularly at the first part at 3.6-6 cm, requires some (small) 
amount of thermal absorption since synchrotron self-absorption (SSA) can never generate a 
spectral index exceeding +2.5 ; and 3) There is no evidence for systematic spectral index 
change of the emission within the available observations. While short term fluctuations 
of the spectral index cannot be ruled out, none appear to occur at the same epoch at all 
frequencies. 



4.4. Brightness temperature 



Figure [3] shows the estimated brightness te mperatures for each frequency versus time 
for SNe 19941 (left) and for 1993J (right; from lweiler et al. 1120071 ). While the precise 
temperature measurements depend on the distance used and the assumed (or, for SN 1993J, 
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measured) shock wave speed, it is clear that the two plots have a systematic difference at 
early times. In the case of SN 19941, the temperatures for each frequency appear to remain 
approximately constant near the maximum possible value of 3 x 10 11 K (see Section 12.31 
until the source becomes optically thin. However, the measured brightness temperatures 
for SN 1993J at the earliest times increase for each frequency until the emission becomes 
optically thin. This is evidence that the dominant absorption mechanism is different 
between these two SNe, with the primary absorption mechanism for SN 1993J not directly 
related to the emission mechanism (i.e., external thermal electrons), while the primary 
absorption mechanism for SN 19941 is physically associated with the emission mechanism 



[i.e., synchrotron electrons). This simple comparison makes it clear 



jhat SSA is dominant 



Weiler et al. 



(120071 ). SN 1993J 



in SN 19941 at early times while, as demonstrated by 
requires significant contributions to the early absorption from both FFA and SSA. Thus, 
SN 19941 is the first observed radio supernova (RSN) to clearly demonstrate the need for a 
SSA dominated component to model its radio light curves. 



4.5. External distant absorption 



The parameter K A (see Section 13 . 1 j) is interpreted as being due to a constant amount 
of ionized material between us and SN 19941 that is sufficiently far from the SN to be 
unaffected by the expanding shock wave. It is relate d to th e emis sion measure (EM) of this 
intervening ionized hydrogen by Equation (1-223) of 



Lang 



(119861 ) so that 



EM = 8.93 x 10 7 K 4 (T c /10 4 K) L35 pc cnT 



(13) 



where T e is the electron temperature of the ionized absorbing medium. Thus, the value for 
K4 in Table H] implies a presence along the line of sight to SN 19941 of ionized hydrogen 
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with an emission measure of EM ~ 4.1 x 10 (T e /10 K) pc cm 



Similar to the case of SN 1978K (IMontes. Weiler. and Panagia 



1997 



Chuetal 



19991) 



this thermal absorption is likely to be due to an intervening H II region and not part of the 
matter from the presupernova mass-loss wind. Most HII regions with an emi ssion measure 



higher than ~ 10 6 pc cm 6 are known to be rather compact sources (see, e.g., 



Wilson et al. 



1970 



Reifenstein et al. 



1970J); for example, the Orion Nebula has an average emission 



19731 ) . If we adopt a similar 



measure of 1.8 x 10 pc cm and a diameter of ~ 2 pc (lAUenl 
diameter of ~ 2 pc for this HII region presumed to be along the line-of-sight to SN 19941, 
we derive an average electron density of 2 x 10 3 pc cm -6 and a total mass of the ionized 
gas of M(HII) ~ 250 M . The volume emission measure is then n 2 e V ~ 5.1 X 10 62 cm -3 , 
which implies a steady source of ionizing radiation of Nl - 1.3 x 10 50 Lyman-continuum 
photo ns s" 1 . Such a flux is higher than that provided even by a bright O-type star (see, 



e.g., 



Panagia 



19731 ) and suggests the presence of a young, compact star cluster in the 



foreground of SN 19941. 

In order not to make the projection of such a powerful cluster in front of SN 19941 too 
geometrically improbable, the cluster should be spatially close to the supernova itself; say, 
separated by no more than ~ 3 pc, which would give a probability of chance projection of 
~ 3%. This then implies that the progenitor of SN 19941 may have been a member of that 
cluster and, by association, was a massive star. 



4.6. Pre-supernova mass-loss rate 

The best-fit value for the parameter K-i (see Section I3.1D is related to the density 
of a presupernova stellar wind. If the density distribution in the circumstellar material 
(CSM) around SN 19941 is proportional to r~ 2 , as appropriate for a constant mass- loss rate, 



constant velocity progenitor wind confirmed by llmmler. Wilson. &: Terashima I (120021 ). our 
measurement of K 2 = 30 in Table H] corresponds to a CSM external thermal absorption as a 
function of time of: 



r(A, t) = 0.35 \- 2A (t/l day)" 3 (14) 

Assuming for the moment that the progenitor star for SN 19941 was a red supergiant 
(RSG), inserting Equation (!T4"|) into equation (16) of Weiler et al. (1986), adopting an 
electron temperature of T e = 20,000 K for the absorbing CSM gas, and assuming a wind 
velocity for this suggested red supergiant progenitor of w = 10 km s _1 , we obtain a mass 
loss rate of: 



M = 1.8 x 10~ 7 {v s /2Q, 000 km s" 1 ) 1 ' 5 (w/10 km s" 1 ) (T e /20, 000 K) 0S8 M yr" 1 (15) 



This estima te of M is much lower than the estim ate of 1.1 x 10 5 (w/10 km s *) M Q yr 1 



obtained by 



Immler. Wilson, fc Terashima 



((2002|) from their analysis of the SN 19941 X-ray 



emission data. 



Clumping of the CSM (or a systematic deviation from spherical symmetry) might 
provide a simple explanation for this apparent discrepancy if most of the ionized gas is 
present in the form of large, isolated structures. In such a case, the mass-loss rate derived 
from the X-ray emission, which is integrated over the whole CSM, would be much higher 
than that estimated from the radio absorption, because the latter arises only from the 



material along the line of sight whic h may not intercept any dense c 



Indeed, in deriving their M estimate 



Immler. Wilson, fc Terashima 



umps in the CSM. 



(120021 ) assumed that 



the CSM was both uniform and spherically symmetric. If such appreciable clumping is 
present, the average density of the emitting material could be grossly underestimated for 
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radio measurements of absorption along the line-of-sight by a factor equal to the square 
root of the clumping factor. The X-ray estimated mass-loss rate (proportional to the ratio 
of the X-ray luminosity to the average density) would then be overestimated by the same 
factor. Conversely, the radio absorption-based estimate would provide only a lower limit 
to the true mass-loss rate. Thus, with a suitable selection of the clumping factor such a 
large difference between the radio and X-ray estimates could be accounted for without 
abandoning the RSG progenitor assumption. For example, if we take a "compromise" value 
of 



\jM(radio) x M(X ray) ~ 1.4 x 1(T 6 (w/10 km s _1 ) M yr" 1 



(16) 



to be appropriate for SN 19941 progenitor's wind, the clumping factor would be ~ 40 and 
the fraction of mass in the diffuse component ~ 15%. 

On the other hand, if the progenitor of SN 19941 was not a RSG star but a Wolf-Rayet 
(WR) star, as has been suggested for Type Ib/c SNe by many authors (see, e.g., 



van den Bergh 



1988 : 



Gal- Yam et al. 



2007 



Smartt 



20091 . and references therei n), the 



Willis 



WR progenitor would have a much higher wind speed (> 10 3 km s -1 ; see, e.g., 
19961 ). as opposed to ~ 10 km s _1 for a RSG star. Thus, a WR progenitor would revise 
the above estimate of the presupernova mass-loss rate from Equation [16] upwards to 
M > 1.4 x 10~ 4 M yr _1 , a m ass-loss rate that is quite possible for a WR star towards the 



end of its evolution (see, e.g., 



Willis 



19961 ). 



5. Conclusions 



We present detailed radio observations of SN 19941 at multiple wavelengths from a 
few days to more than 8 years after explosion. Although there have been many radio 
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observations of Type Ib/c supernovae since SN 19941, SN 19941 remains the best studied 
and most heavily sampled example of the class. Also, even though observations since those 
of SN 19941 have indicated that synchrotron self-absorption can play a role in the early 
radio light curves of supernovae that have relatively low density circumstellar matter, such 
as Type Ib/c and Type lib, SN 19941 was the first example that was so well sampled 
early enough that such a conclusion was inescapable. Also, as seen in other RSNe (see 
Section I4.2[) this low density circumstellar thermal material in SN 19941 may be irregularly 
distributed as discussed in Section I4T61 

As with previous detailed observations of the radio emission from supernovae, there 
is no sign of systematic spectral index evolution with time. However, a changing spectral 
index with frequency, rather than with time, is a distinct possibility, with the highest 
observed frequencies perhaps indicating a somewhat flatter spectral index. 

As has been noted in the few other radio supernovae which have been sufficiently well 
studied, the radio light curves, even though their general properties evolve smoothly and are 
well described by existing models, show significant "bumps and wiggles." SN 19941 shows a 
significant local "dip" between ~ 10 to ~ 55 days after explosion that is interpreted as a 
temporary decrease in the average density of the circumstellar medium, perhaps implying 
the origin of SN 19941 in an exploding Wolf-Rayet star in a wide binary system with an 
early-type stellar companion, although changing mass-loss rate from the progenitor star 
cannot be ruled out. 

The usual scenario for core-collapse (Type Ib/c and Type II) supernovae is that the 
progenitor was a red supergiant. However, based on the information available from these 
new observations, it appears that a Wolf-Rayet star undergoing a period of very high 
mass-loss rate is the most likely progenitor star for SN 19941. 
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Table 1. Centimeter Flux Density Measurements for SN 19941 



Obs. 


Age 


Tel. 


S(20 cm) ± rms 


S(6 cm) ± rms 


S(3.6 cm) ± rms 


S(2 cm 


) ± rms 


S(1.3 cm) ± rms 


Date 


(days) 


Config. 


(mjy) 


(mjy) 


(mjy) 


(m 


Jy) 


(mjy) 


31-Mar-94 


= 0.00 


















03-Apr-94 


3.167 


VLA-A 






0.554 ± 0.062 








2.580 ± 0.388 


04-Apr-94 


4.120 


VLA-A 






0.737 ± 0.120 


2.780 ± 0.469 


3.707 ± 0.891 


05-Apr-94 


5.350 


VLA-A 


<0.090 a 


0.477 ± 0.091 


1.743 ± 0.108 


6.300 ± 0.522 


8.299 ± 0.891 


06-Apr-94 


6.040 


VLA-A 




0.700 ± 0.115 












06-Apr-94 


6.135 


Camb 








10.800 


± 


1.472 




06-Apr-94 


6.483 


VLA-A 




1.400 ± 0.115 












07-Apr-94 


7.028 


VLA-A 






2.412 ± 0.152 










07-Apr-94 


7.070 


VLA-A 






3.295 ± 0.190 










07-Apr-94 


7.190 


Camb 








13.000 


± 


2.385 




07-Apr-94 


7.462 


VLA-A 






3.656 ± 0.194 


13.100 


± 


1.062 


16.300 ± 1.676 


08-Apr-94 


8.020 


Camb 








14.500 


± 


2.086 




08-Apr-94 


8.288 


VLA-A 


<0.207 a 


1.250 ± 0.153 


4.990 ± 0.250 


14.698 


± 


1.202 


15.458 ± 1.673 


09-Apr-94 


9.180 


Camb 








16.000 


± 


2.561 




09-Apr-94 


9.386 


VLA-A 




2.020 ± 0.144 


6.514 ± 0.334 


17.200 


± 


1.409 


17.100 ± 1.810 


10-Apr-94 


10.561 


VLA-A 


<0.225 a 


2.123 ± 0.176 


8.930 ± 0.455 


19.171 


± 


1.488 


18.593 ± 2.028 


ll-Apr-94 


11.180 


Camb 








17.500 


± 


2.658 




12-Apr-94 


12.180 


Camb 








17.200 


± 


2.282 




12-Apr-94 


12.395 


VLA-A 


<0.132 a 


4.120 ± 0.221 


11.43 ± 0.577 


15.860 


± 


1.249 


12.400 ± 1.349 


13-Apr-94 


13.220 


Camb 








13.900 


± 


2.436 




13-Apr-94 


13.372 


VLA-A 


<0.162 a 


4.643 ± 0.271 


12.605 ± 0.646 


17.178 


± 


1.365 


12.610 ± 1.414 


14-Apr-94 


14.220 


Camb 








15.600 


± 


2.536 




14-Apr-94 


14.367 


VLA-A 


0.300 ± 0.067 


5.053 ± 0.267 


13.000 ± 0.660 


16.500 


± 


1.382 




15-Apr-94 


15.220 


Camb 








15.000 


± 


2.500 




15-Apr-94 


15.360 


VLA-A 


<0.258 a 


5.587 ± 0.292 


14.500 ± 0.733 


16.500 


± 


1.364 


10.500 ± 1.178 


16-Apr-94 


16.220 


Camb 








16.500 


± 


2.593 




16-Apr-94 


16.356 


VLA-A 


0.370 ± 0.095 


6.451 ± 0.338 


16.200 ± 0.820 


15.600 


± 


1.236 


11.000 ± 1.252 


17-Apr-94 


17.210 


Camb 








16.100 


± 


2.200 




17-Apr-94 


17.363 


VLA-A 


0.463 ± 0.105 


7.664 ± 0.400 


17.200 ± 0.871 


16.900 


± 


1.328 


10.600 ± 1.185 


18-Apr-94 


18.199 


VLA-A 




7.390 ± 0.392 


15.770 ± 0.810 


15.743 


± 


1.215 


10.590 ± 1.205 


18-Apr-94 


18.210 


Camb 








15.900 


± 


2.555 





Table 1 — Continued 



Obs. Age Tel. S(20 cm) ± rms S(6 cm) ± rms S(3.6 cm) ± rms S(2 cm) ± rms S(1.3 cm) ± rms 

Date (days) Config. (mJy) ( m Jy) ( m Jy) ( m Jy) ( m Jy) 



19-Apr-94 


19.220 


Camb 








14.000 ± 2.441 




20-Apr-94 


20.200 


Camb 








13.500 ± 2.018 




20-Apr-94 


20.494 


VLA-A 


0.421 ± 0.107 


9.250 ± 0.493 


16.483 ± 0.850 


8.757 ± 0.960 


5.673 ± 1.185 


21-Apr-94 


21.200 


Camb 








14.000 ± 2.052 




22-Apr-94 


22.210 


Camb 








13.000 ± 1.985 




23-Apr-94 


23.056 


VLA-A 


0.386 ± 0.089 










23-Apr-94 


23.150 


Camb 








12.000 ± 1.921 




24-Apr-94 


24.336 


VLA-A 


0.688 ± 0.125 


11.130 ± 0.568 








25-Apr-94 


25.140 


Camb 








10.000 ± 1.414 




25-Apr-94 


25.333 


VLA-A 


0.720 ± 0.122 


10.760 ± 0.551 


12.300 ± 0.651 




4.620 ± 0.816 


28-Apr-94 


28.163 


VLA-A 


0.683 ± 0.115 


13.713 ± 0.693 


14.723 ± 0.760 


7.503 ± 0.733 


4.990 ± 0.706 


29-Apr-94 


29.200 


Camb 








10.000 ± 1.281 




30-Apr-94 


30.320 


VLA-A 


1.618 ± 0.167 


14.063 ± 0.726 


14.650 ± 0.757 


8.720 ± 0.848 


6.293 ± 0.971 


02-May-94 


32.140 


Camb 








9.800 ± 1.400 




03-May-94 


33.130 


Camb 








9.100 ± 1.212 




03-May-94 


33.279 


VLA-BnA 


1.170 ± 0.142 


16.823 ± 0.847 


11.440 ± 0.591 




3.988 ± 0.624 


10-May-94 


40.228 


VLA-BnA 


2.660 ± 0.267 


17.550 ± 0.880 


13.040 ± 0.656 


5.330 ± 0.477 


2.180 ± 0.555 


12-May-94 


42.090 


Camb 








8.000 ± 1.281 




15-May-94 


45.179 


VLA-BnA 


2.470 ± 0.270 


16.995 ± 0.861 


12.788 ± 0.652 


6.555 ± 0.615 


3.533 ± 0.662 


16-May-94 


46.100 


Camb 








8.000 ± 1.281 




17-May-94 


47.100 


Camb 








8.400 ± 1.306 




19-May-94 


49.252 


VLA-BnA 


3.420 ± 0.343 


15.940 ± 0.798 


11.040 ± 0.554 


5.610 ± 0.523 


2.810 ± 0.539 


23-May-94 


53.208 


VLA-BnA 




15.860 ± 0.809 


11.305 ± 0.570 


5.380 ± 0.547 


2.680 ± 0.722 


24-May-94 


54.060 


Camb 








5.500 ± 1.141 




27-May-94 


57.060 


Camb 








5.500 ± 0.890 




28-May-94 


58.136 


VLA-BnA 


4.660 ± 0.466 


14.030 ± 0.704 


8.400 ± 0.425 


4.240 ± 0.473 


2.820 ± 0.645 


Ol-Jun-94 


62.188 


VLA-BnA 


5.080 ± 0.508 


12.400 ± 0.626 


7.670 ± 0.392 


3.390 ± 0.441 


<2.250 a 


17-Jun-94 


78.165 


VLA-B 


7.190 ± 0.729 


9.640 ± 0.497 


5.860 ± 0.310 


2.920 ± 0.404 


1.910 ± 0.356 


29-Jun-94 


90.173 


VLA-B 


10.940 ± 1.103 


7.620 ± 0.388 


4.240 ± 0.221 


2.420 ± 0.359 


2.050 ± 0.596 


07-Jul-94 


98.027 


VLA-B 


10.530 ± 1.065 


7.340 ± 0.369 


4.000 ± 0.207 


2.150 ± 0.323 


<1.620 a 


08-Aug-94 


130.937 


VLA-B 


9.410 ± 0.951 


5.020 ± 0.262 


2.750 ± 0.150 


1.480 ± 0.348 


<2.400 a 
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Table 1 — Continued 



Obs. 


Age 


Tel. 


S(20 cm) ± rms 


S(6 cm) ± rms 


S(3.6 cm) ± rms 


S(2 cm) ± rms 


S(1.3 cm) ± r 


Date 


(days) 


Config. 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


22-Aug-94 


144.959 


VLA-B 


11.480 ± 1.155 


3.760 ± 0.234 


1.930 ± 0.154 


1.620 ± 0.476 


<1.980 a 


Ol-Scp-94 


154.753 


VLA-B 


9.250 ± 1.277 


2.073 ± 0.174 


1.415 ± 0.139 


<1.320 a 


<4.050 a 


30-Scp-94 


183.774 


VLA-CnB 


10.510±1.196 


2.824±0.213 


1.530 ± 0.097 


<1.830 a 


<2.160 a 


13-Oct-94 


196.687 


VLA-C 


<2.841 a 


2.203±0.264 


1.117±0.169 


<1.140 a 


<1.917 a 


07-Nov-94 


221.513 


VLA-C 




1.257 ± 0.191 


1.220 ± 0.117 


<0.990 a 


<2.970 a 


05-Dec-94 


249.708 


VLA-D 




2.209 ± 0.357 


1.024 ± 0.176 


<2.220 a 




05-Jan-95 


280.507 


VLA-D 








<1.650 a 




06-Apr-95 


371.173 


VLA-D 








<1.032 a 


<1.920 a 


16-Jun-95 


442.147 


VLA-DnA 


2.397 ± 0.253 


0.855 ± 0.102 








06-Oct-95 


554.790 


VLA-BnA 


1.087 ± 0.121 


0.503 ± 0.113 


0.391 ± 0.072 


<0.636 a 




12-Jan-96 


652.500 


VLA-CnB 




0.460 ± 0.046 




<0.480 a 




05-Oct-96 


919.847 


VLA-DnA 


0.580 ± 0.152 


0.250 ± 0.071 


<0.390 a 






13-Jun-99 


1900.000 


VLA-DnA 


<0.222 a 










05-Apr-02 


2927.000 


VLA-A 


0.150 ± 0.026 











a All upper limits are 3<r 
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Table 2. Millimeter Flux Density Measurements for SN 19941 



Obs. 
Date 


Age 
(days) 


Tel. 


Flux Density (S) ± rms* 
(mJy) 


Frequency 
(GHz) 




Reference 




31-Mar-94 


= 0.00 














05-Apr-94 


5.45 


OVRO 


11.0 ± 0.6 


99 


Phillips, 


Chandler, & Rupen 


(1994) 


06-Apr-94 


6.45 


OVRO 


<12 b 


218 


Phillips. 


Chandler, & Rupen 


(1994) 



a Quoted errors are taken from the original references. 
b All upper limits are 3ct 
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Table 3. VLA Calibrator Flux Density Measurements' 1 



Obs. S(20 cm) S(6 cm) S(3.6 cm) S(2 cm) S(1.3 cm) 



Date 


(mjy) 


(mjy) 


(mjy) 


(mJy) 


(mjy) 


03-Apr-94 






0.539 b 




0.338 b 


04-Apr-94 






1.014 c 


0.953 c 


0.830 c 


05-Apr-94 


0.904 


1.038 


1.014 


0.953 


0.830 


06-Apr-94 




1.020 c 








07-Apr-94 






0.994 


0.959 c 


0.880 c 


08-Apr-94 


0.899 c 


1.025 c 


1.000 c 


0.954 c 


0.875 c 


09-Apr-94 




1.040 c 


1.019 c 


0.980 c 


0.910 c 


10-Apr-94 


0.894 


1.013 


0.987 


0.956 


0.920 


12-Apr-94 


0.887 c 


1.006 c 


0.974 c 


0.940 c 


0.864 c 


13-Apr-94 


0.884 


1.017 


0.993 


0.955 


0.865 


14-Apr-94 


0.875 c 


1.026 


1.027 c 


0.976 c 




15-Apr-94 


0.875 c 


1.017 c 


1.029 c 


0.993 c 


0.839 c 


16-Apr-94 


0.875 c 


1.017 c 


1.029 c 


0.993 c 


0.839 c 


17-Apr-94 


0.875 c 


1.060 c 


1.032 c 


0.976 c 


0.810 c 


18-Apr-94 




0.984 


0.963 


0.930 


0.901 


20-Apr-94 


0.875 c 


1.005 c 


0.964 c 


0.918 c 


0.773 


23-Apr-94 


0.875 c 










24-Apr-94 


0.875 c 


1.015 c 








25-Apr-94 


0.904 


1.015 


0.946 




0.689 


28-Apr-94 


0.866 


0.995 


0.937 


0.881 


0.683 


30-Apr-94 


0.887 c 


0.993 c 


0.958 c 


0.890 c 


0.792 c 


03-May-94 


0.905 


0.992 


0.940 




0.900 


10-May-94 


0.890 


0.946 


0.930 


0.859 


0.846 


15-May-94 


0.899 


0.983 


0.949 


0.900 


0.816 


19-May-94 


0.868 


0.948 


0.920 


0.830 


0.764 


23-May-94 




0.986 


0.933 


0.863 


0.797 


28-May-94 


0.872 


0.976 


0.939 


0.878 


0.822 


Ol-Jun-94 


0.894 


0.996 


0.963 


0.874 


0.905 


17-Jun-94 


0.894 


0.991 


0.957 


0.875 


0.793 


29-Jun-94 


0.883 


0.958 


0.918 


0.829 


0.786 


07-Jul-94 


0.882 


0.927 


0.885 


0.783 


0.753 
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Table 3 — Continued 







eff: rrn \ 




0/9 nrr) \ 




Date 


(mjy) 


(mjyj 


(mjyj 


(mjyj 


(mjy) 


08-Aug-94 


0.881 


0.912 


0.879 


0.846 


0.827 


22-Aug-94 


0.926 


0.863 


0.784 


0.635 


0.633 


Ol-Sep-94 


0.906 c 


0.852 c 


0.807 c 


0.657 


0.665 


30-Sep-94 


0.885 


0.831 


0.830 


0.654 


0.528 


13-Oct-94 


0.840 c 


0.834 c 


0.822 c 


0.705 c 


0.610 c 


07-Nov-94 




0.826 


0.812 


0.755 


0.693 


05-Dcc-94 




0.832 


0.782 


0.680 




05-Jan-95 








0.726 




06-Apr-95 








0.697 


0.634 


16-Jun-95 


0.760 


0.764 








06-Oct-95 


0.705 


0.750 


0.738 


0.702 




12-Jan-96 




0.726 




0.697 




05-Oct-96 


0.622 


0.775 


0.792 






13-Jun-99 


0.498 










05-Apr-02 


0.537 d 











a Unlcss otherwise noted, the flux density measurements are for J1418+546. 
b Flux density measurements are for J1415+463 
c Interpolated from nearby measurements 
d Flux density measurement is for 1327+434 



-36 - 



Table 4. Best fit parameters for SN 19941 and Comparison with SN 1993J. 



Parameter 


SN 19941 


SN 1993J a 


Ki 


4.78 x 10 3 


4.8 x 10 3 


a 


-1.22 


-0.81 


P 


-1.42 


-0.73 


K 2 


3.0 x 10 1 


1.6 x 10 2 


5 


= -3.00 b 


-1.88 


K 3 




4.3 x 10 5 


8' 


— 


-2.83 


K 4 


4.55 x 1CT 2 


=0 


K 5 


6.90 x 10 5 


2.62 x 10 3 


8" 


-4.08 


-2.05 


K 6 


9.96 x 10 3 


=0 


8"' 


-2.63 




Distance (Mpc) 


8.9 


3.63 


Date of explosion (to) 


1994 March 31 


1993 March 28 


Time to S 6 cm pcak (days) 


34.4 


133 


Se cm peak ( m Jy) 


16.2 


96.9 


Lq cm pcak (erg s _1 Hz -1 ) 


1.5 x 10 27 


1.5 x 10 27 


M (M yr" 1 ) 


> 1.4 x 10~ 4c 


0.5-5.9 x 10~ 6d 



a Weiler et al. 2007 

b Adopted. Because the amount of thermal material surrounding 
SN 19941 is so low, it is not possible to determine 8 and K3 from 
the data fitting and nominal values were adopted. 

c Under the assumption of a Wolf-Rayet progenitor; see Section 
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d Weiler et al. (2007) found that the mass-loss rate from 
SN f993J increased from ~ 5.4x 10~ 7 Af yr _1 at the time of shock 
breakout to - 5.9 x 10~ 6 M Q yr^around day - 3100 (~ 8000 yr 
before explosion). 
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Fig. 1. — The radio light curves for SN 19941 for: top left 0.3 cm, top right 1.2 cm, middle 
left 2 cm, middle right 3.6 cm, lower left 6 cm, lower right 20 cm. The lines on all figures 
represent the best fit model as described in the text with the parameters listed in Table 4. 

Fig. 2. — The spectral index evolution for SN 19941 for: top left 1.2-2 cm, top right 2-3.6 
cm, lower left 3.6-6 cm, lower right 6-20 cm. The lines on all figures represent the best fit 
model as described in the text with the parameters listed in Table 4. 

Fig. 3. — (left) The brightness temperature (Tb) evolution, uncorrected for thermal absorp- 
tion, for SN 19941 for, top left 1.2 cm, top right 2 cm, middle left 3.6 cm, middle right 6 cm, 
and lower left 20 cm and (right) t he brightness tempe rature (Tb) evolution, uncorrected for 



thermal absorption, for SN 1993J ( jWeiler et al. 



20071 ) for, from left to right, 0.3 cm (cross) 



1.2 cm (filled square), 2 cm (open square), 3.6 cm (filled circle), 6 cm (open triangle), 20 cm 
(filled triangle), 49 cm (star), and 90 cm (open diamon d). The horizontal dashed lines in 



all plots denote the limiting value of Tb — 3 x 10 K (iKellermann Sz Pauliny-Toth 



1969 



Readhead 



199J). 



Fig. 4. — Cartoon (taken from lWeiler et al. I (120021 )). not to scale, of the supernova and its 
shocks, along with the stellar wind established circumstellar medium (CSM), the interstellar 
medium (ISM), and more distant ionized hydrogen (HII) absorbing gas. The radio emission 
is thought to arise near the blastwave front. The expected locations of the several absorbing 
terms in Equations ([5] [12]) are illustrated. 

Fig. 5. — (left) Log of the deviations of the ratio of observed to model flux density for 
SN 19941 for the best fit model listed in Table H] from 1.2 cm (top) to 20 cm (bottom) and 
(right) composite of the deviations at all wavelengths shown on the left, with a smoothed 
curve added to guide the eye. 
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